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§ ; ABSTRACT 



Previously published packings of equal disks in an equilateral triangle have dealt with up 
to 21 disks. We use a new discrete-event simulation algorithm to produce packings for up 
to 34 disks. For each n in the range 22 < n < 34 we present what we believe to be the 



r^ ' densest possible packing of n equal disks in an equilateral triangle. For these n we also list the 

a ■ 

second, often the third and sometimes the fourth best packings among those that we found. 

In each case, the structure of the packing implies that the minimum distance d{n) between 
►^ . disk centers is the root of polynomial P„ with integer coefficients. In most cases we do not 

ly-! I explicitly compute Pn but in all cases we do compute and report d{n) to 15 significant decimal 

\Q ■ digits. 

,—1. I Disk packings in equilateral triangles differ from those in squares or circles in that for 



o 
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in this case. Based on our computations, we present conjectured optimal packings for seven 
other infinite classes of n, namely 



triangles there are an infinite number of values of n for which the exact value of d(n) is known, 
namely, when n is of the form A(/c) := -^ — -. It has also been conjectured that d{n — l) = d{n) 



X 

H 
. .C^ ; n = A(2A;) + 1, A(2/c + l) + l,A(/c + 2)-2, A(2A: + 3)-3, 

A(3/c + 1) + 2, 4A(A;), and 2A(A; + 1) + 2A(A;) - 1 . 

We also report the best packings we found for other values of n in these forms which are larger 
than 34, namely, n = 37, 40, 42, 43, 46, 49, 56, 57, 60, 63, 67, 71, 79, 84, 92, 93, 106, 112, 121, 
and 254, and also for n = 58, 95, 108, 175, 255, 256, 258, and 260. We say that an infinite 
class of packings of n disks, n = n(l), n(2), ...n(k), ..., is tight , if [l/d{n{k) + 1) — l/d(n{k))] is 
bounded away from zero as k goes to infinity. We conjecture that some of our infinite classes 
are tight, others are not tight, and that there are infinitely many tight classes. 



1 Introduction 

Geometrical packing problems have a long and distinguished history in combinatorial math- 
ematics. In particular, such problems are often surprisingly difficult. In this note, we describe 
a series of computer experiments designed to produce dense packings of equal nonoverlapping 
disks in an equilateral triangle. It was first shown by Oler in 1961 [O] that the densest packing 
of n = A{k) := —^ — - equal disks is the appropriate triangular subset of the regular hexagonal 
packing of the disks (well known to pool players in the case of n = 15). It has also been con- 
jectured by Newman [N] (among others) that the optimal packing of A(fc) — 1 disks is always 
obtained by removing a single disk from the best packing for A(fc), although this statement 
has not yet been proved. The only other values of n (not equal to A(fc)) for which optimal 
packings are known are n = 2, 4, 5, 7, 8, 9, 11 and 12 (see Melissen [Ml], [M2] for a survey). 

As the number n of packed disks increases, it becomes not only more difficult to prove 
optimality of a packing but even to conjecture what the optimal packing might be. In this 
paper, we present a number of conjectured optimal packings. These packings are produced on 
a computer using a so-called "billiards" simulation algorithm. A detailed description of the 
philosophy, implementation and applications of this event-driven algorithm can be found in [L], 
[LS]. Essentially, the algorithm simulates a system of n perfectly elastic disks. In the absence 
of gravitation and friction, the disks move along straight lines, colliding with each other and the 
region walls according to the standard laws of mechanics, all the time maintaining a condition 
of no overlap. To form a packing, the disks are uniformly allowed to gradually increase in size, 
until no significant growth can occur. Not infrequently, it can happen at this point that there 
are disks which can still move, e.g., disk 3 in t7al3 (see Fig. 1.1). 

Every packing of n disks occurring in the literature for n different from A(/c) and A{k) — 1 
which has been conjectured or proved to be optimal was also found by our algorithm. These 
occur for n = 13, 16, 17, 18, and 19 (see [Ml], [MS]). This increases our confidence that the 
new packings we obtain are also optimal. The new packings cover two "triangular periods": 
21 = A(6) to A(7) to A(8) = 36. 

In addition, we conjecture optimal packings for seven infinite classes of n, namely, n = 
A{2k)+1, A{2k+1)+1, A{k+2)-2, A{2k+3)-3, A(3/c+l)+2, 4A(A;), and 2A{k+l)+2A{k)-l, 
where k = 1, 2.... Each class has its individual pattern of the optimal packings which is different 
from patterns for other classes. These were suggested by the preceding packings, and we give 




t7al3 

0.366025403784439 




13 bonds 



t7al6 

0.366025403784439 16 bonds 



Figure 1.1: Two equivalent but nonisomorphic densest packings of 7 disks. 

packings for some additional values of these forms, namely, n = 37, 40, 42, 43, 46, 49, 56, 57, 
60, 63, 67, 71, 79, 84, 92, 93, 106, 112, 121, and 254, as well as for n = 58, 95, 108, 175, 255, 
256, 258, and 260. 

We say that an infinite class of packings of n disks, n = n{l),n(2), ...n{k), ..., is tight , if 
[l/d{n(k) + 1) — l/d(n(k))] is bounded away from zero as k goes to infinity. We conjecture 
that some of our infinite classes are tight, others are not, and that there are infinitely many 
tight classes. 

2 The packings 

We performed a small number of runs with n = 21, 27, 28, 35 and 36 disks. In every case, 
the resulting packings were consistent with the existing results (n = A{k)) and conjectured 
(n = A(/c) — 1). The bulk of our efforts concentrated on the other 11 values of n, for 21 < n < 
36. These are presented in Figures 3.1 to 3.11. 

To navigate among the various packings presented we will use the labeling system illustrated 
by Fig. 3.1 t22a. Here, n = 22, "a" denotes that the packing is the best we found, "b" would 
be the second best (as in t23b in Fig. 3.2), "c" would be third best, and "d" would be fourth 
best. 



Small black dots in the packing diagrams are "bonds" whose number is also entered by 
each packing. For example, there are 47 bonds in t22a. A bond between two disks or between 
a disk and a boundary indicates that the distance between them is zero. The absence of a 
bond in a spot where disk-disk or disk-wall are apparently touching each other means that the 
corresponding distance is strictly positive, though perhaps too small for the resolution of the 
drawing to be visible. For example, there is no bond between disk 1 and the left side of the 
triangle in tl8a (Fig. 2.2); according to our computations, the distance between disk 1 and the 
side is 0.0048728... of the disk diameter. (Packing tl8a was constructed in [Ml].) Each disk in 
most of the packings is provided with a label which uniquely identifies the disk in the packing. 
This labeling is nonessential; it is assigned in order to facilitate referencing. 




tl7a40 

0.211324865405187 40 bonds 



tl7a42 

0.211324865405187 42 bonds 



tl7a43 

0.211324865405187 43 bonds 




tl7b36 tl7b42ns tl7b42s 

0.208735129275750 36 bonds 0.208735129275750 42 bonds 0.208735129275750 42 bonds 



Figure 2.1: The best (tl7a40, tl7a42, tl7a43) and the next-best (tl7b36, tl7b42ns, tl7b42s) 
packings of 17 disks. 



Each disk normally has at least three bonds attached. The polygon formed by these bonds 
as vertices contains the center of the disk strictly inside. This is a necessary condition for 
packing "rigidity". In [LS], where the packing algorithm was applied to a similar problem, the 
disks without bonds were called "rattlers." A rattler can move freely within the confines of 
the "cage" formed by its rigid neighbors and/or boundaries. (If we "shake" the packing, the 
rattler will "rattle" while hitting its cage.) t22a has two rattlers, disks 3 and 5. In the packing 
diagrams, all disks, except for the rattlers, are shaded. 

A number with 15 significant digits is indicated for each packing in the figures, e.g., the 
number 0.17939 69086 11866 for packing t22a. This number is the disk diameter d{n) which is 
measured in units equal to the side of the smallest equilateral triangle that contains the centers 
of all disks. For packing t22a such a triangle is the one with vertices at the centers of disks 22, 
17, and 12. This unit of measure for d{n) conforms with previously published conventions. 

Sometimes several packings exist for the same disk diameter. An example is t7al3 and 
t7al6 in Fig. 1.1. Thus, we distinguish such packings by suffixing their labels with the number 
of bonds. Other examples are tl7a40, tl7a42, and tl7a43 in Fig. 2.1, t22b42 and t22b50 in 
Fig. 3.1. However, even the number of bonds may not distinguish different packings of the 
same disk diameter; for example, tl7b42ns and tl7b42s in Fig. 2.1, where the provisional "ns" 
stands for "non-symmetric" and "s" for "symmetric." 

We point out that the a-packings of 17 and 18 disks that we show have previously been 
given by Melissen and Schuur [MS], who also conjecture their optimality. 

3 Additional comments 

Fig. 3.2: Two more c-packings for 23 disks that are not shown in the figure were generated: 
t23c55.1 and t23c55.2. Both have 55 bonds. t23c55.1 can be obtained by combining the left 
side of t23c53 with the right side of t23c57. t23c55.2 is a variant of t23c55.1. 

Fig. 3.3: Disk 20 in t24c56 and in t24c59 is locked in place because its center is strictly 
inside the triangle formed by the three bonds of disk 20. In both packings, the distance of the 
disk center to the boundary of this enclosing triangle is the distance to the line between bonds 
with the left side of the triangle and disk 24, and is 0.0317185... of the disk diameter. 

Fig. 3.4: The given d-packing of 25 disks t25d60 is symmetric with respect to the vertical 
axis. An equivalent non-symmetric d-packing t25d53 was also obtained in which all disks are 
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located in the same places as in t25d60, except for disks 5, 12, 13, 14, 23, and 24. These six 
disks form a pattern which is roughly equivalent to that formed by disks 10, 14, 19, 25, 20, 
and 22, respectively, in t25b. Disk 24 in t25d53 is a rattler. 

Fig. 3.6: Only one of the two b-packings of 29 disks we found is shown, namely, t29b63.2. 
The other b-packing, t29b63.1, differs in the placements of only disks 2, 3, 4, 7 as explained in 
Section 4. 

Fig. 3.8: Four a-packings of 31 disks exist; only three are shown in the figure; the fourth 
one, t31a81.1, is described in Section 4. 

Fig. 3.10: In t33a, the gap between disk 8 and left side is 0.0017032... of the disk diameter. 
In t33c, disk 7 is stably locked by its bonds with 3, 6, and 29. However, the distance from disk 
7 center to the line on bonds with disks 3 and 6 is only 0.0002097575.... of the disk diameter. 
As a result, the cage of rattler disk 5 in t33c is very tight: the gap between disk 22 and disk 5 
or disk 18 and disk 5 does not exceed 4 x 10^^ of the disk diameter. 

Fig. 3.11: In t34a, the small gaps between "almost" touching pairs disk-disk or disk- wall 
take on only three values (relative to the disk size): in pairs 20-31, 16-26, 23-27, 18-19, 1-27 
the gap is 0.021359..., in pairs left-32, right-29 it is 0.024750..., and in pairs 4-34, 7-22, it is 
0.042561... Similarly, there are only three values of gaps in each of t34b, t34c, and t34d. 

t34b: in pairs 18-19, 23-27, 17-28, 20-31, 16-26 the gap is 0.019583...; in pairs 

left-32, right-29 it is 0.022686...; in pairs 4-34, 7-22, it is 0.039035... 

t34c: in pairs 12-17, 22-27, 3-10, 14-21, 4-34. 3-16 the gap is 0.018864...; in pair 

left-15 it is 0.021850...; in pair 19-24 it is 0.037606... 

t34d: in pairs 2-4, 26-32, 15-22, 12-21, 3-16, 7-16 the gap is 0.018681...; in pair 

left-27 it is 0.021637...; in pairs 13-33, 19-30 it is 0.037242... 




tl8a 

0.203465240539124 40 bonds 




tl8b36 

0.203464834591373 36 bonds 




tl8b40 

0.203464834591373 40 bonds 




tl8b43 

0.203464834591373 43 bonds 



Figure 2.2: The best (tl8a) and the next best (tl8b36, tl8b40, tl8b43) packings of 18 disks. 




t22a 

0.179396908611866 47 bonds 




t22b42 

0.179132453213560 42 bonds 




t22b50 

0.179132453213560 50 bonds 




t22c 

0.178763669382058 46 bonds 



Figure 3.1: The best (t22a), the next-best (t22b42, t22b50), and the third-best (t22c) packings 
of 22 disks. 




t23a 

0.175153309170525 52 bonds 




t23b 

0.174962364462008 52 bonds 




t23c53 




t23c57 



0.174457630187009 53 bonds 



0.174457630187009 57 bonds 



Figure 3.2: The best (t23a), the next-best (t23b), and the third-best (t23c53, t23c57) packings 
of 23 disks. 
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t24a 

0.174457630187010 63 bonds 




t24b 

0.171024411616889 53 bonds 




t24c56 

0.170613243353863 56 bonds 




t24c59 

0.170613243353863 59 bonds 



Figure 3.3: The best (t24a), the next-best (t24b), and the third-best (t24c56, t24c59) packings 
of 24 disks. 
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t25a 

0.169065874417891 63 bonds 




t25b 

0.167753380810744 58 bonds 




t25c 

0.167714758299772 57 bonds 




t25d60 



0.167685455574367 60 bonds 



Figure 3.4: The best (t25a), the next-best (t25b), the third-best (t25c), and a fourth-best 
(t25d60) packings of 25 disks. 
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t26a 

0.166738399395271 63 bonds 




t26b 

0.166732017260692 63 bonds 




t26c 

0.166698218018074 63 bonds 




t26d 

0.166695724000574 62 bonds 



Figure 3.5: The best (t26a), the next-best (t26b), the third-best (t26c), and the fourth-best 
(t26d) packings of 26 disks. 
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t29a 




t29b63.2 



0.152189614060732 63 bonds 



0.152172645377571 63 bonds 




t29c 

0.152161553910934 65 bonds 




t29d 

0.152109020552728 63 bonds 



Figure 3.6: The best (t29a), a next-best (t29b63.2), the third-best (t29c), and the fourth-best 
(t29d) packings of 29 disks. 
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t30a 

0.150761500215428 73 bonds 




t30b 

0.149057883638389 71 bonds 




t30c 

0.149047199036657 68 bonds 




tSOd 

0.148951506290246 67 bonds 



Figure 3.7: The best (t30a), the next-best (tSOb), the third-best (tSOc), and the fourth-best 
(tSOd) packings of 30 disks. 
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t31a79 




0.148543145110506 79 bonds 



t31a81.2 

0.148543145110506 81 bonds 




t31a82 




t31b 



0.148543145110506 82 bonds 



0.146207330881888 71 bonds 



Figure 3.8: The best (t31a79, t31a81.2, t31a82) and the next-best (t31b) packings of 31 disks. 
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t32a 

0.145102169183849 76 bonds 




t32b 

0. 1449847274688 12 72 bonds 




t32c 




t32d 



0. 1 4465 1911 950450 7 1 bonds 



0.144616419018845 77 bonds 



Figure 3.9: The best (t32a), the next-best (t32b), the third-best (t32c), and the fourth-best 
(t32d) packings of 32 disks. 
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t33a 




t33b75 



0.143447408371201 79 bonds 



0.143447385418276 75 bonds 




t33b79 

0.143447385418276 79 bonds 




t33c 

0.143309997215537 74 bonds 



Figure 3.10: The best (t33a), the next-best (t33b75, t33b79), and the third-best (t33c) packings 
of 33 disks. 



18 




t34a 

0.142869646754496 84 bonds 




t34b 

0. 14286764768 1 844 83 bonds 




t34c 

0. 14286688784583 1 84 bonds 




t34d 

0.142866698669904 83 bonds 



Figure 3.11: The best (t34a), the next-best (t34b), the third-best (t34c), and the fourth-best 
(t34d82) packings of 34 disks. 



19 

4 Conjectures for individual packings 

Each packing diagram we give can imply several different statements: 

(I) There exists a valid configuration of nonoverlapping disks with all pairwise distances 
marked by bonds equal to zero, and those not marked by bonds strictly positive, and 
with disk diameter equal to the indicated value with a relative error of less than 10^^^. 

(II) The configuration is rigid: no disk or set of disks except for rattlers can be continuously 
displaced from the indicated positions without overlaps. 

(Ill) The configurations are correctly ranked. That is, the a-packing really is optimal, the 
b-packing is second best, etc. 

We believe (I) and (II) are correct. As to (III), we hope the statement is correct with 
respect to the a-packings. In other words, we believe these are the optimal packings. We are 
less confident for the lower ranked packings. For example, if someone finds a new packing 
in between our c- and d-packings, we will not be astounded. We provide these mainly for 
comparison purposes, and to serve as benchmarks for other packing algorithms. 

It would also not be surprising to discover a nonisomorphic packing to one we have presented 
which has exactly the same disk diameter and the same number of bonds (e.g., as in tl7b42ns 
and tl7b42s in Fig. 2.1). 

5 Conjectures for infinite classes 

Dense packings in an equilateral triangle seem to "prefer" to form blocks of dense triangles 
and arrangements that are nearly so. In this section we describe seven infinite classes where 
we think we have found the optimal packings. Each class has its individual pattern of the 
optimal packings, which is different from the patterns for other classes. However, since they 
are the result of the particular packings we found, which themselves are only conjectured to 
be optimal, then the general conjectures have even less reliability. We still think they might 
serve as useful organizers for the maze of published dense packings. 

4A(k). The best packing (we found) of 24 = 4A(3) disks in Fig. 3.3 consists of four triangles, 
each with A(3) = 6 disks. The best packings of 12 disks (in [Ml]) in Fig. 5.1, and even 4 disks in 
Fig. 4.1 have the same form. The packings we obtained while experimenting with 40 = 4A(4) 
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t4a 

0.577350269189626 9 bonds 



t8a 

0.343070330817254 18 bonds 



Figure 4.1: The best packing of 4 disks (t4a) and 8 disks (t8a). 

and 60 = 4A(5) disks (see Fig. 5.2), and also with 84 = 4A(6) and 112 = 4A(7) disks have 
the same structure as well. 

A simple analysis of the patterns obtained implies that d{4A{k)) 



2k-2+V3' 

If we fit members of class 4A(.) within the boundaries of the triangular periods, i.e., among 
members of the class A(.), then every other period has exactly one n of the form 4A(A;) lying 
almost exactly at the middle of the period. 



2A(k + 1) + 2A(k) — 1. For each k there are k + 1 distinct best packings: two for 7 disks 
(Fig. 1.1), three for 17 disks (Fig. 2.1), four for 31 disks (three of these four are shown in 
Fig. 3.8), and five for 49 disks (four of these five are shown in Fig. 5.3). 
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tlla tl2a tl3a 

0.275255128608411 22 bonds 0.267949192431123 30 bonds 0.251813236653061 30 bonds 

Figure 5.1: The best packings of 11 disks (tlla), of 12 disks (tl2a), and of 13 disks (tl3a). 





t40a 

0.129331793710034 108 bonds 



t60a 

0.102753265449690 154 bonds 



Figure 5.2: The best packings of 40 disks (t40a) and 60 disks (t60a). 
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t49al30 

0.114520634618068 130 bonds 




t49al32.2 

0.114520634618068 132 bonds 




t49al33 




0.114520634618068 133 bonds 



t49al32.3 

0.114520634618068 132 bonds 



Figure 5.3: Four (out of the five existing) best packings of 49 disks. 
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t37a 

0.132739196122383 84 bonds 




t56a 

0.105229325597427 133 bonds 



Figure 5.4: The best packings of 37 disks (t37a) and 56 disks (t56a). 

There are three different numbers of bonds in these packings; the smallest number of bonds 
is in the packing with a rattler, then k — 2 packings each of which has a "cavity" and the same 
number of bonds, and finally one more packing without a rattler or a cavity with the largest 
number of bonds. The four triangles, two small and two large, that illustrate the expression 
2A{k + 1) + 2A{k) - 1, can be seen in the packings with a rattler (tl7a40, t31a79, t49al30). 
The two larger triangles are defective: both coalesce a corner disk (disk 6 in tl7a40, disk 4 in 
t31a79, disk 4 in t49al30). Packing t31a81.2 is obtained from t31a79 by the left larger triangle 
acquiring its corner disk 4 and pushing disks 31 and 29 from the left side of the other large 
triangle down into the cavity formed . If we push only disk 31, we obtain a packing t31a81.1 
(not shown). If we push all three disks 31, 29, 10 into the cavity (and rotate the resulting 
structure to recover the symmetry with respect to the vertical axis), we obtain the fourth best 
packing t31a82. 

A simple analysis of the patterns obtained implies that d{2A{k+l) + 2A{k) — l) = ^r— — y=. 
Our experiments with 71 = 2A(6) + 2A(5) — 1 disks produced the same patterns for the best 
packing. 

If we fit the members of class 2A{k + 1) + 2A(/c) — 1 into the boundaries of the triangular 
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tl6a33.1 



tl6a33.2 



tl6b 



0.216227269309782 33 bonds 0.216227269309782 33 bonds 0.215843184970192 32 bonds 

Figure 5.5: The best (tl6a33.1, tl6a33.2) and the next-best (tl6b) packings of 16 disks. 

periods, as we did for the class 4A{k), we find that every other period has exactly one n of 
the form 2A{k + 1) + 2A(fc) — 1 which lies almost exactly in the center of the period. Thus, 
classes 2A(A;-|-1)-|-2A(A;) — 1 and 4A(fc) are "parity-complementary" to each other. Beginning 
with the second triangular period 3 to 6, each period has exactly one term of one or the other 
class; even periods contain terms of the class 4A(fc) and odd periods contain terms of the class 
2A{k + 1) + 2A{k) -1. 

The pattern of one of the parity-complementary class pair can be obtained from the pattern 
of the other class by a simple transformation. For example if we eliminate the eight bottom 
disks 2, 42, 6, 25, 15, 29, 32, 47, and the rattler 4 in t49al30 (Fig. 5.3), we obtain the pattern 
of t40a (Fig. 5.2). 



A(2k) -|- 1. As we noted earlier, A{m) densely packed disks in an equilateral triangle form 
a perfect hexagonal lattice with m disks on a side. When m = 2A; is even, the structure 
with A(2fc) -|- 1 disks adjusts itself to one extra disk as follows. The top 2k — 1 rows remain 
packed hexagonally, and the bottom row ripples to accommodate 2A; -|- 1 disks instead of 2k. 
In this ripple of the bottom row, the 1st, 3rd, ... {2k + l)th disk beginning from the left corner 
remain attached to the bottom, while the 2nd, 4th, ... 2A;th disk rise and attach themselves 
to 2nd, 4th, ... 2A;th disks respectively, of the row above. Thus, k — 1 rigid cages are formed. 
The k — 1 disks of the row above to which no disk is attached from below fall off into these 
cages and become rattlers. The first seven terms of the class A(2fc) -|- 1 are: t4a (k = 1, no 
rattlers since A; — 1 = 0, see Fig. 4.1), tlla (constructed in [Ml] with one rattler; see Fig. 5.1), 
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t46a 

0. 1 17209943988839 106 bonds 




t46bl06.2 

0. 1 17208402974392 106 bonds 



Figure 5.6: The best (t46a) and a next-best (t46bl06.2) packings of 46 disks. 

t22a (see Fig. 3.1) with two rattlers, t37 and t56 (Fig 5.4) with 3 and 4 rattlers, respectively, 
t79a (194 bonds, 5 rattlers, d(79) = 0.0871159038791759), and tl06a (267 bonds, 6 rattlers, 
(i(106) = 0.0742982999063026). We do not reproduce the diagrams here for the latter two 
packings; their patterns are identical to the class description given above. 

A(2k + 1) + 1. When m = 2A; + 1, A; = 1, 2..., the odd parity of m causes a more complex 
adjustment to the extra disk. The bottom row ripples in a non-symmetric way; the ripple 
creates k cages for rattlers and a cavity; see t29a (Fig. 3.6) for k = 3. 

Notice that packing t29b63.2 (Fig. 3.6) has almost the same structure as t29a, except for 
the cage that consists of disks 2, 3, 7, 9, 6, 5, and 8, is depressed and disk 4 in t29b63.2 is not a 
rattler, and a nonrattler 6 in t29a becomes a rattler in t29b63.2. The same two modifications 
exist for k = 4 (i.e., n = 46), and the modification with the depressed cage, t46bl06.2, is 
again inferior (Fig. 5.6). Beginning with k = 5 (i.e., n = 67), while both modifications exist, 
they exchange their roles: the depressed one becomes the best, t67al61.2, while the other one 
becomes the inferior one, t67b (Fig. 5.7). For example, t92a228.2 is the modification with 
the depressed cage, while t92b is the other one (Fig. 5.8). The same pattern is displayed by 
packing tl21a307.2 (307 bonds, 6 rattlers, d(121) = 0.0691630188894699), for which we omit 
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t67al61.2 




0.0952213722318656 161 bonds 



t67b 

0.0952208239815187 161 bonds 



Figure 5.7: A best (t67al61.2) and the next-best (t67b) packings of 67 disks. 

the diagram here. 

Labels t29b63.2, t46bl06.2, t67al61.2, t92a228.2, and tl21a307.2 have the suffix 2 in 
them because there exist equivalent packings t29b63.1, t46bl06.1, t67al61.1, t92a228.1, and 
tl21a307.1, respectively. The latter differ from the former in the placement of only 4 disks. 
An easy way to explain this is to look at the second term of the class n = A{2k + 1) + 1 for 
n = 16 disks (Fig. 5.5). In this case both modifications exist and both deliver the optimum, 
tl6.a33.1 and tl6.a33.2. They differ in the placement of disks 2, 3, 4, and 7. 

The side rattler disk 15 in tl6b can be considered a precursor for the side rattler 29 in 
t29d. The same side-rattler pattern was observed in lower ranked packings for k = 4 (n = 46), 
k = 5 {n = 67), k = 6 {n = 92), and k = 7 {n = 121). 

Classes A{2k) + 1 and A(2/c + 1) -|- 1 are a parity complementary pair, similar to the pair 
of classes 4A(A;) and 2A{k -|- 1) -|- 2A(/c) — 1 considered above. 



A(k -|- 2) — 2. While the optimal packings of A{k -|- 2) — 1 disks are always (apparently) 
perfectly hexagonal with a single disk removed, the removal of two disks from a hexagonal 
arrangement is never optimal. The first two terms t4a and t8a (Fig. 4.1) suggest no common 
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t92a228.2 



0.0801405249366249 228 bonds 



t92b 

0.0801398309712360 228 bonds 



Figure 5.8: A best (t92a228.2) and the next-best (t92b) packings of 92 disks. 

pattern. Looking at the next case tl3a (Fig. 5.1) suggests the pattern for fc > 3 of a packed 
triangle A(A;) at the top, supported by two sparse rows of disks, each of which lacks a disk 
compared to what would be there in a perfect hexagonal packing. The top-A(/c)-plus-two- 
sparse-rows packing indeed exists for any k > 3 and is rigid. In particular, the pattern appears 
again in the best packings for k = 4 (tl9a in Fig. 5.11). 

However, the optimality of this pattern does not continue for /c > 4, as can be seen in 
Fig. 3.5 where the best packing t26a has a different pattern. The top-A(/i;)-plus-two-sparse- 
rows pattern is not even among the top four packings for n = 26. The pattern for t26a and 
t26b persists for the next term (see t34a and t34b in Fig. 3.11) but then roles become reversed 
for 43 disks (see t43a and t43b in Fig. 5.10). Will the pattern of packing t43a remain optimal 
for larger values of k? Unfortunately, our algorithm fails to obtain stable packings for 53 (or 
larger values of A (A; + 2) — 2) disks. 



A(3k + 1) + 2 =A(3k - 1) + (2k + 1) A(2). 30 = A(5) + 5A(2) and packing t30a (Fig. 3.7) 
can be viewed as a A (5) triangle on top from which disk 8 fell off and became a rattler, 
supported by five triangles A (2) from below; t30a is the second term of the class. To produce 
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t42a 

0.125451671189752 108 bonds 



t63a 

0.100080300212798 165 bonds 



Figure 5.9: The best packings of 42 disks (t42a) and of 63 disks (t63a). 

this structure for the next value of k we add two triangles A (2) on the bottom and three more 
layers of disks to enlarge the top triangle to become a A (8) (and, in general, this procedure 
is repeated for larger values of k). Indeed our experiments with n = 57 (A; = 3) and n = 93 
{k = 4) did produce this structure in the best packings (see t57a and t93a in Fig. 5.12). For 
A: = 1 we have n = 12, a degenerate case with no rattlers. The number of rattlers in this 
packing for general k is k — 1. 

A(2k + 3) - 3 =A(2k) + (2k + 1) A(2). The first term is packing tl2a (Fig. 5.1) which also 
belongs to class 4A(.). The second term is packing t25a (Fig. 3.4) which can be viewed as 
a A(2A;) triangle on top supported by 2/c + 1 alternating A(2) triangles below. This pattern 
of the second term, whose description also fits the first term, is more apparent in the third 
and fourth terms (see t42a and t63a shown in Fig. 5.9). According to our experiments the 
next terms t88a and tll7a do not continue this pattern but give way to patterns which are 
somewhat similar to those of the class A (A; + 2) — 2 considered above. (We omit diagrams for 
both t88a and tll7a.) 
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6 How good are the packings? 

Let us compare the values of the best packings with the only bound currently available, 
namely one based on the inequality of Oler [O] . This inequality has the following form (see [FG] 
for a simple proof): Let K he a, compact, convex subset of E^ with area A(K) and perimeter 
P{K). li p{K) denotes the maximum number of points that can be placed in K so that any 
pair has mutual distance at least 1, then the following inequality holds: 

(1) p{K) < -^AiK) + ^PiK) + 1 . 

Inverting (1) and applying it with p{K) = n and K being an equilateral triangle of side length 
L{n), we obtain 

(2) L{n) > -{-3 + VS^^+T) := t{n) . 

In Fig 6.1 we plot the difference 

5{n) := L{n) - t{n) 

versus n for selected values of n < 121. A dot with a circle around it indicates that the 

corresponding value has been proved to be optimal; a dot without a surrounding circle or an 

open square indicates that the value is only conjectured to be optimal. For up to 37 disks 

there were only four values which did not fall into one of our infinite classes, namely, n = 18, 

23, 32, and 33; those are indicated by open squares. The values which have been associated 

with classes are connected by lines, with a distinct type of line for each class. 

We should point out that for each n, the value of the largest disk diameter d{n) and the 

value of L{n) are reciprocally related, i.e., d{n)L{n) = 1. Thus, L{-^ — -) = k — 1 for k > 1. 

If our conjecture for n = A.{k) + 1 = -^ — - + 1 is correct then it would follow that 

/k(k + l) \ 2 

lim 6 — + 1 < ^ - 1 = 0.1547 .... 

fc-.oo V 2 / ~ ^3 

Using the explicit values for d{4A{k)) and d{2A{k + 1) + 2A(/c) — 1) given in Section 5, 
we have L(4A(/c)) = 2A; - 2 + ^/3 and L{2A{k + 1) + 2A{k) - 1) = 2k - 1 + ^/S, from 
which it follows that both lim^-^oo (5(4A(/c)) and limfc_^oo S{2A{k + 1) + 2A(/c) — 1) are at most 
-v/3 — 3/2 = 0.2321 . . . with the distance between fcth term and the limit being of the order 
of l/k. 

We conjecture that for each of the classes n = A{2k) + 1, A(2A; + 1) + 1, and A{3k + 1) + 2 
the value of 6{n) is bounded away from zero. In fact, we believe that for any fixed c > 0, 
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the value of 6{A{k) + c) is bounded away from zero. In other words, packings of A(A;) disks 
are so tight that any attempt to accommodate even one additional disk noticeably worsens 
the packing in that L{n) increases by at least some positive amount independent of n. In this 
sense the packings for the class A(A;) are "tight". 

On the other hand, we believe that after any fixed positive number of disks are added to 
A(fc) disks, any other fixed number of disks can be added without substantial "damage" to 
6{n) (asymptotically). Thus, for example, it would seem that if limfe_»oo S{n{k)) exists for each 
of the class n{k) = A{2k) + 1, A{2k + !) + !, and A{3k + 1) + 2 (and the limits probably do 
exist), then all three limits are equal. In this sense the classes A{2k) + 1, A(2fc + 1) + 1, and 
A(3fc + l) + 2 are "loose". 

Similarly, the classes n = A{k + 2) — 2 and A(2fc + 3) — 3 are "loose" in the sense that 
we can add one disk to the best (conjectured) packing without noticeable change of 6{n) for 
sufficiently large n. This follows from the fact that limfc_^oo (5(A(/i;) — p) ^ for p fixed. The 
latter limit is obvious by noticing that a lower-bounding packing for A{k) — p disks when k is 
sufficiently large is simply the densest packing of A(/c) disks with p disks removed; for such a 
packing, 5 is asymptotically 0. 

Formally, we say that an infinite class of packings of n disks, n = n(l),n(2), ...n(k), ..., is 
loose, if lim.k^^[5{n(k) + 1) — 5{n{k))] = 0. Because we believe this limit exists for any class 
we consider, each class has to be either tight or loose. 

We further conjecture that the classes 4A(A;) and 2A(A; + 1) + 2A(A;) — 1 are tight, similar 
to the class A{k). 

Are there other tight classes? Here is our argument in favor of the existence of a countable 
infinity of distinct tight classes. We believe that if each densest packing of n disks for n = 
n{l),n{2),...n{k),.., consists of a fixed number, say r, of densely packed triangles A(.), then 
6{n{k) + 1) — 6{n{k)) is bounded away from zero as k goes to infinity. 

Thus, we have the following task: for each r from some infinite set find a sequence 
n{l),n{2),...n{k),.., so that the densest packing of n{k) disks for all k has the "same pat- 
tern" and consists of exactly r densely packed triangles. Note that a priori we are not able to 
define what the "same pattern" is (and hence we have no formal definition of what a "class" 
is); but after producing a class the pattern is usually clear. 

Let us consider a two-parameter family of numbers n = np{k), p = 1,2..., k = 1,2..., of 
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the form 

(3) np(fc) = A((A; + l)(p + 1) - 2) + A; = A((A; + l)p - 1) + {2p + 1)A(A;) . 

Two equal expressions for np{k) are given in (3). The second expression suggests r = 2{p + 1) 
triangles. If we take, for example, k = p = 2, we get ^2(2) = 30 disks and r = 6. The 
conjectured t30a (Fig 3.7) indeed consists of 6 densely packed triangles, if we attach rattler 8 
to the top triangle A (5). Take now 58 = n2(3). The pattern of our experimental packing t58a 
(Fig 6.2) looks like t30a (Fig 3.7) with 6 triangles again (with the rattler attached to the top 
triangle). 

Thus, just as t30a is a member of the class A(3A; — 1) + {2k + 1)A(2) for k = 2, t58a is 
(perhaps) a member of the class A(4A; — 1) + {2k + 1)A(3) for k = 2. The pattern of the kth. 
packing of this class is composed of 2{k + 1) densely packed triangles. Is this class tight or 
loose? We believe it is loose because A(4/c — 1) + {2k + 1)A(3) = A(3/c + 1) + 2, and a class of 
the form "A(.) + const" is always loose (we think). Incidentally, the number of the triangles 
in the class is unbounded with k. 

However, if (as we believe) the sequence of densest packings np{k),k = 1,2,3,... for fixed 
p = 2 can be continued with all packings having the same pattern of six triangles, then t58a 
might also be a member of the class A(2A; + 1) + 5A(A;) for A; = 3. The next term in the latter 
class would be the densest packing of n = n2(4) = A(9) + 5 A (4) = 95 disks. Our experiments 
with 95 disks, indeed, produced the desired pattern of six triangles in the densest packing t95a 
(Fig. 6.2). This reinforces our suspicion that the class n2{k), k = 1,2, ..., exists in which each 
densest packing consists of six densely packed triangles. The class n2{k), k = 1,2, ..., should 
be tight because each densest packing in it consists of a fixed number of triangles. 

By increasing p we are moving into a different class, which is again tight if the conjecture 
above is correct. Thus, for p = 3 we have the sequence 71-3(1) = 22,113(2) = 57, 71,3(3) = 
108, 113(4) = 175, .... Packings t22a (Fig 3.1), t57a (Fig 5.12) indeed each consist of 2(p+l) = 8 
densely packed triangles. Our experiments with 108 and 175 disks yield the same pattern in 
the best packings (Fig 6.2) so the class n^{h) probably exists too. In the same way, the class 
np{k), k = 1,2, ... exists for any fixed index p and has a distinct pattern with r = 2{p + 1) 
triangles and p — 1 rattlers. 

If this is correct, then Figure 6.2 can be seen as the 2x2 submatrix for 2 < p < 3 and 
3 < /c < 4 of the matrix of dense packings of np{k) disks where 1 < k,p < 00. By traversing a 
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row or a column of this matrix we obtain a distinct infinite class of packings. Our conjecture 
is that each row class is tight and each column class is loose. 

This matrix contains three infinite classes conjectured in Section 5: the row at p = 1 is the 
class 4A(fc), the column at fc = f is the class A(2p) + 1, and the column at /c = 2 is the class 
A{3p + 1) + 2 = A{3p - 1) + (2p + 1)A(2). 

The conjecture about the full matrix is also reinforced by the fact that k and p in (3) are 
unique for each given value of n = np{k). This can be easily seen using the first expression for 
np{k) in (3). To further test our matrix conjecture, we generated the list of all n of the form 
n = np{k) for n < 300. These are: 4, 11, 12, 22, 24, 30, 37, 40, 56, 57, 58, 60, 79, 84, 93, 95, 
106, 108, 112, 137, 138, 141, 144, 172, 174, 175, 180, 192, 196, 211, 220, 254, 255, 256, 258, 
260, 264, and 280. Some increments in this increasing sequence are small. Specifically, in each 
following subsequence increments do not exceed 2: (11, 12), (22, 24), (56, 57, 58, 60), (93, 95), 
(106, 108), (137, 138), (172, 174, 175), and (254, 255, 256, 258, 260). 

Now, take for example, 255 = nj{2) and 256 = ^5(3). These are two "almost" equal 
numbers of disks. However according to the matrix conjecture they should produce different 
patterns of densest packings: the pattern for 255 should consist of one large and 15 small 
triangles with 6 rattlers and the pattern for 256 of one large and 11 small triangles and 4 
rattlers. Similarly, the matrix conjecture prescribes specific patterns for the densest packing of 
the other numbers of disks n of this sequence, e.g., 254 = nii(l), 258 = ^3(5), and 260 = ^2(7). 

One might think it would be a stress test for both the matrix conjecture and our packing 
procedure to try to pack these numbers of disks. Note that many of the packings for smaller 
values in the sequences above have been generated (as discussed above) and they all conform 
to the matrix conjecture. Thus, we experimented with packing n = 254, 255, 256, 258, and 
260 disks. Recall that the procedure of packing has no idea, so to say, of the desirable packing. 
Starting with random initial conditions the disks perform chaotic movements, they collide with 
each other and with the boundaries millions of times and each collision evaluation is subject 
to roundoff error. 

The experiments turned out to be not so difficult. (Case of 53 disks proved to be harder.) 
As expected, the best packing of 254 disks has the pattern of its class /S.[2k) + 1 with 10 rattlers 
((i(254) = 0.0467170396481042, 679 bonds; we omit the diagram). Fig.6.3 shows the patterns 
of packings t255a, t256a, 1258, and t260a. These too are consistent with formula (3). Note 
that because of the large number of disks in the packings the scale of drawing in Fig.6.3 is 
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small and bonds are not seen. The pictures with a larger scale (omitted here) show that all 
the bonds exist in right places. 
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Figure 6.1: Discrepancy between side length of the triangle and its lower bound for different 
n. 
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t58a 

0.103809894100183 151 bonds 




t95a 

0.0791578042660103 256 bonds 




tl08a 



0.0736689924964931 287 bonds 



tl75a 



0.0569014905693285 479 bonds 



Figure 6.2: The best packings of 58 disks (t58a), 95 disks (t95a), 108 disks (tlOSa), and 175 
disks (tl75a). 
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t255a 

0.0465999333293566 693 bonds 



t256a 

0.0465050212813994 703 bonds 





t258a 

0.0463515918933559 719 bonds 



t260a 

0.0462256883510349 733 bonds 



Figure 6.3: The best packings of 255 disks (t255a), 256 disks (t256a), 258 disks (t258a), and 
260 disks (t260a). 
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t43a 

0. 125004280307977 108 bonds 




t43b 

0.125003232541680 108 bonds 




t43c 

0. 125002529873623 108 bonds 




t43d 

0. 12500223482957 1 106 bonds 



Figure 5.10: The best (t43a), the next-best (t43b), the third-best (t43c) and the fourth-best 
(t43d) packings of 43 disks. 
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tl9a 



tl9b 



tl9c 



0.200321458983439 45 bonds 0.200200894290590 45 bonds 0.200186574270937 45 bonds 



Figure 5.11: The best (tl9a), the next best (tl9b), and the third best (tl9c) packings of 19 
disks. 





t57a 

0.104447019187802 143 bonds 



t93a 

0.07978572868437232 240 bonds 



Figure 5.12: The best packings of 57 disks (t57a) and 93 disks (t93a) 
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7 Discussion 

While a finite number of patterns for infinite classes have been tentatively identified to date 
(two one-parameter patterns known or conjectured previously joined by several such patterns 
in Sec. 5 and a two-parameter "matrix" pattern in Sec. 6) a countable infinity of such patterns 
and classes probably exists. Furthermore, each value of n may well be a member of one or 
more such classes. Thus, the values n = 18, 23, 32, and 33, which were not placed into classes 
in this paper, may well be members of as yet unidentified classes of packings with complex 
patterns. In fact, a fixed value of n may be on the paths of many, possibly infinitely many, 
such classes. 12 disks gives an example of this: it is on the path of the class 4A{k) and it is 
also the first term of the classes A(3A; + 1) + 2 and A(2fc + 3) — 3. As the value of n increases 
along the path of a class, "hesitations" of the best pattern may occur, wherein several different 
nonequivalent patterns coexist among the rigid packings and compete for the title of the best. 
A resolved case of such hesitation occurs for the class A(2fc + 1) + 1 where for k > 5 (n > 67) 
two equivalent best patterns finally emerge (at least according to our experiments). We were 
not able to confirm by experiments the winning pattern for the class A(fc + 2) — 2. Will such 
hesitation always be resolved in favor of one of the competing patterns in a finite initial segment 
of the path? 
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